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Abstract 

DNA is an ideal candidate to organize matter on the nanoscale, primarily due to the specificity 



^4— I ' and complexity of DNA based interactions. Recent advances in this direction include the self- 

O ! 
. ' assembly of colloidal crystals using DNA grafted particles. In this article we theoretically study 



the self-assembly of DNA-caged particles. These nanoblocks combine DNA grafted particles with 



T^ [ more complicated purely DNA based constructs. Geometrically the nanoblock is a sphere (DNA 

Q ' grafted particle) inscribed inside a polyhedron (DNA cage). The faces of the DNA cage are 



open, and the edges are made from double stranded DNA. The cage vertices are modified DNA 
junctions. We calculate the equilibriuim yield of self-assembled, tetrahedrally caged particles, and 



(N 
> 

r~^ I discuss their stability with respect to alternative structures. The experimental feasability of the 

a^' 

^f) • method is discussed. To conclude we indicate the usefulness of DNA-caged particles as nanoblocks 



in a hierarchical self-assembly strategy. 



I. INTRODUCTION 

DNA is one of the most celebrated tools in the nanoscience toolbox. This approach 
was pioneered in the laboratory of N. Seeman, where some of the first schemes for build- 
ing nanostructures from specially engineered oligonucleotide sequences were proposed. A 
number of objects have been successfully constructed, including DNA cubes, multiple armed 
DNA junctions, DNA crystals, and DNA knots (_|l| , [2] , |3| , [^ ) . There have been several re- 
cent experimental advances in this direction ( 5|,[6|), including the encapsulation of a single 



molecule inside a DNA tetrahedron 



u. 



There has also been a surge of interest in utilizing 
the specific interactions of complementary single-stranded DNA (ssDNA) to organize par- 
ticles on the nanoscale. One recent advance in this direction is the self-assembly of three 



dimensional body centered cubic crystals from DNA grafted nanoparticles 



Up until 



this point, most of the studies reported the formation of small clusters or random aggregation 
of particles, as opposed to the self-assembly of ordered structures (|lO|],[ll[, 12] )■ 

The potential complexity of DNA based interactions provides a means to design signifi- 
cantly more complicated nanoblocks. In this paper we propose a method to self-assemble 
DNA-caged particles (see Fig. [T]). These nanoblocks are composite materials which are con- 
structed by combining DNA grafted nanoparticles with specially designed DNA sequences. 
Geometrically the nanoblock is a sphere (DNA grafted nanoparticle) inscribed inside of a 
polyhedron (DNA cage). The polyhedron faces are open, and the cage edges are made of 
double-stranded DNA (dsDNA). The cage vertices are modified DNA junctions ([l3|, |l4|). 
Each vertex of the cage carries a unique ssDNA sequence available for hybridization. This 
vertex "coloring" makes these nanoblocks ideal candidates as building blocks for hierarchical 
self-assembly strategies. 

The plan for the paper is the following. We first introduce of our self-assembly proposal. 
Details of the proposal are discussed for a particular implementation in which the DNA cage 
is a regular tetrahedron. We theoretically calculate the melting profile for the DNA cage 
self-assembled around the DNA grafted particle. We demonstrate an equilibrium regime in 
which the DNA-caged particle is the dominant structure in solution, and discuss its stability 
with respect to alternative structures. We conclude by discussing how DNA-caged particles 
could be used as the building blocks in a hierarchical self-assembly strategy. 
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FIG. 1: (Color online). Graphical depiction of the scheme for self-assembling DNA-caged particles. 
The cage edges are constructed from dsDNA rods with terminal ssDNA sequences on either end. 

II. SELF-ASSEMBLY PROPOSAL 



In this section we discuss the details of our self-assembly proposal. The proposal is to 
self-assemble DNA-caged particles by combining DNA grafted nanoparticles with rod-like 
DNA linkers. There is only one type of particle, i.e. all of the ssDNA grafted onto the 
nanoparticle surface have the same nucleotide sequence. The rod-like DNA linkers are 
dsDNA, but each end of the rod terminates in a ssDNA sequence. The dsDNA rods can 
bend significantly when their length L is comporable to the persistence length Ip ~ 50nm 
for dsDNA [15i]. In what follows we consider the case L <^ Ip and treat the dsDNA as 
rigid rods. There are n types of rod-like DNA linkers, since the terminal ssDNA sequences 
on each linker are unique. The particle and linkers are all of the components necessary 
for the self-assembly proposal. The number of types of linkers is determined by the cage 
architecture, in general there will be one type for each edge of the cage. We now turn to 
discuss how the cage can be assembled from the DNA linkers. 

For the sake of concreteness we will consider a particular implementation of this idea in 




FIG. 2: (Color online). The vertex architecture is a modified 4 arm DNA junction. The arrow- 
heads label the 3' end of the ssDNA arms. Note that the portion of arm 1 which is normally 
complementary to arm 4 is missing. Hence the 3' end of arm 4 provides the vertex with a unique 
"color", i.e. a ssDNA sequence available for hybridization. 



which the DNA cage is a regular tetrahedron. In this case there are n = 6 types of DNA 
rods, one for each edge of the tetrahedron. 

These rods can be joined to assemble the cage in the following manner. To construct each 
vertex of the tetrahedron, four ssDNA sequences must be joined. Three of these ssDNA 
sequences are the terminal ssDNA sequences of the rod-like DNA linkers. The fourth ssDNA 
sequence comes from the ssDNA grafted onto the particle surface, which binds the particle 
to the cage. The DNA architecture that accomplishes this task is known as a four arm DNA 
junction (see Fig. [2]). These junctions have been studied extensively, and the sequences can 
be designed so that the vertex is stable 13|, uA, ll6| . 

The problem is now to assign particular sequences to the terminal ssDNA sequences of 
the DNA rods which result in the desired tetrahedral cage, taking into account the proposed 
vertex architecture. One such assignment is proposed below in Fig. [31 We now provide an 
explanation of how Figs. [2] and [3] can be read together to understand the cage construction. 




FIG. 3: (Color online). The proposed sequence assignment to the ssDNA ends of the DNA rods 
which self-assembles into the tetrahedral cage. If S denotes a dsDNA spacer, the sequence assign- 
ments for the 6 rods are: Ri = A2 — S — B2, R2 = B4 — S — C^, R3 = C2 — S — D2, Ri = A^ — S — Di, 
i?5 = B3 — S — Di, Rq = A3 — S — C3. The detailed structure of the vertex architecture (the 
circles in this diagram) is presented in Fig. [5J 



Examine vertex A in Fig. [31 We can see that rods Ri, R4, and Rq are joined together 
at this vertex. Let An denote the nucleotide sequence which plays the role of arm n (see 
Fig. [2]) in vertex A, where n G {1, 2, 3, 4}. The sequence assignments in the caption of Fig. 
[3] tell which rod provides each arm of the DNA junction. For example, Ri = A2 — S — B2 
means that rod -Ri provides arm 2 of vertex A and arm 2 of vertex B. For vertex A, rod -Ri 
provides arm 2, R4 provides arm 4, and Rq provides arm 3. Only arm 1 remains, which is 
provided by the ssDNA grafted onto the particle. In addition, since all of the ssDNA grafted 
onto the particle have the same sequence, by performing this enumeration procedure for each 
vertex we can see that the following four sequences are identical: Ai = Bi = Ci = D^. 

Since it may be difficult (e.g. for steric reasons) to introduce the particle into the fully 
assembled cage, we would like for the particle to assist in the cage building process. This 
has been explicitly taken into account in the sequence designation process. Note that rod 



monomers: 
dimers: / 

trimers: 
tetramers: 

5-mers: 

6-mers: 



y\ 



\ 




4 



\ 








FIG. 4: (Color online). The subsets of the DNA cage (tetrahedron) which can be formed in the 
absence of the particle. For each topologically distinct diagram only one variety is shown. 

i?4 cannot bind at vertex D in the absence of the particle, since it hybridizes to arm 3 of 
the vertex (which conies from the particle). 

With the basic framework in hand, the next task is to determine the relative abundance 
of the various structures that form in a solution of DNA linkers and DNA-grafted nanopar- 
ticles. A similar type of analysis has been performed in our related work on DNA-grafted 
nanoparticles ( 17|], ISJ]). In the next section we calculate the equilibrium yield for a variety 
of these structures. If the self-assembly process is experimentally feasible we should be able 
to demonstrate a regime in which our nanoblock, a single particle surrounded by a fully 
assembled DNA cage, is the dominant structure in solution. 



III. DNA CAGE MELTING 

We first determine the melting profile for DNA cages in the absence of the nanoparticle. 
By taking the proposed vertex numbering scheme (see Fig. [3]) into account, we can enumerate 
all of the possible DNA structures which can form in solution (see Fig. Hj). 



Let Cj denote the concentration of linker Ri. In what follows Co = IM is the standard 
reference concentration. In addition, we use natural units where the Boltzmann constant 
ks = 1. In equilibrium the chemical potential of the various phases will be the same. 
The chemical potential has a contribution from the entropy of dilution, and the effective 
hybridization free energy for creating the the DNA connections at the vertex. For example, 
consider the reaction in which DNA linker rods -Ri and Rq hybridize to form a dimer. 
Equilibrating the chemical potentials yields the following equation: 

Here eig is the free energy for the formation of the dimer pair from the two monomers. In 
this case Eiq is simply the hybridization free energy associated with joining rods Ri and Rq 
together at vertex A. The resulting concentration of the dimer cig is then 

ci6 = ^exp(^^ (2) 

The total concentration of dimers C2 is determined by considering all of the possible dimer 
varieties. 

C2 = C16 + C46 + Ci5 + C25 + C26 + C36 + C35 (3) 

The dimers can be considered freely-jointed rigid rods, owing to the flexibility of the DNA 
junctions which constitute the vertex. 

In general the free energy Sij is equal to the hybridization free energy AGij = AHij — 
Tlog{ASij) for joining rods Ri and Rj. These free energies will depend on the particular 
choice of the DNA nucleotide sequences {Ai,Bj,Ck,Di}. In what follows e denotes the 
average hybridization free energy {sij). 

The reasoning for the trimer structures is largely the same. To write the hybridization 
free energies for the n-mers compactly, we label them by the indices for the rods which 
constitute the structure. For example, for the trimer composed of rods Ri, R2, and R^ the 
effective free energy is written £125. Looking at Fig. [31 £125 can be decomposed into a sum 
of hybridization free energies for joining two arms at a vertex, i.e. £125 = £15 + £25- The 
same decomposition can be done for all of the n-mers. We adopt the same notation for the 
concentration of the structures. The concentration of the trimer C125 formed by the reaction 
R1 + R2 + R5 is 

C1C2C5 I —£l25\ /,>, 

ci25 = -^^ exp y-^^j (4) 



For some of the DNA structures there is one additional comphcation. For any diagram 
which contains a closed loop, we must calculate the change in configurational entropy as- 
sociated with making the connection which closes the loop. In these cases we calculate 
the overlap density Ce// which relates the effective hybridization free energy e to the bare 



hybridization free energy e ( 18| , 19| ) . 



e.p|_£j.^exp(Z£| (5) 

^ /P(ri,r)P(r2,r)rf-V 
""'^^ (/P(r,rOrf3r)' ^ ^ 

Here P(r, r') is the probability distribution for the chain of DNA linkers which starts at r' 
and ends at r. The canonical example is the conversion of a trimer which is a chain of three 
freely-jointed links into a closed triangle. For rigid DNA linkers each of length L the result 
is quite simple. 

'^^f = 8^ (^) 

Details for the calculation are provided in an appendix. 

Continuing the enumeration procedure for the tetramers, 5-mers, and 6-mer, we can 
write down expressions for the concentration of all the DNA structures which can form in 
the absence of the particle. Writing down the equations for conservation of DNA linkers 
results in a system of 6 coupled polynomial equations of order 6 in the concentrations of 



monomers Cj. Here c*°* is the total initial concentration of linkers of type j. 

cf* = Ci + Ci5 + Ci6 + Ci25 + Ci46 + Ci26 + C135 + C136 + C156 + C1235 + Cl236 (8) 

+Cl246 + Ci346 + C1256 + C1356 + C1456 + C12346 + Cl2356 + Ci2456 + C13456 + C123456 
C2 = C2 + C25 + C26 + C125 + C236 + Ci26 + C235 + C246 + C256 + C1235 + C1236 + C1246 (9) 

+C2346 + Ci256 + C2356 + C2456 + Ci2346 + Ci2356 + Ci2456 + C23456 + Ci23456 
4°* = C3 + C35 + C36 + C236 + C135 + Cl36 + C235 + C346 + C356 + Cl235 + Cl236 + Cl346 (10) 

+C2346 + C2356 + C1356 + C3456 + C12346 + Ci2356 + C13456 + C23456 + C123456 
C4 = C4 + C46 + Ci46 + C246 + C346 + Ci246 + C1346 + C2346 + C1456 + C2456 + C3456 (11) 

+C12346 + Ci2456 + C13456 + C23456 + C123456 
4°* = C5 + Ci5 + C25 + C35 + C125 + C135 + Ci56 + C235 + C256 + C356 + Cl235 + Cl256 (12) 

+C2356 + Cl356 + C1456 + C2456 + C3456 + Cl2356 + Ci2456 + C13456 + C23456 + Ci23456 
Cq = Cg + C16 + C26 + C36 + C46 + Ci46 + C236 + C126 + C136 + C156 + C246 + C256 + C346 (13) 

+C356 + Ci236 + Ci246 + C1346 + C2346 + Ci256 + C2356 + C1356 + C1456 + C2456 

+C3456 + C12346 + C12356 + Ci2456 + C13456 + C23456 + Ci23456 

By solving these equations for the monomer concentrations we can plot the melting profile 
(see Fig. [5]). The plot is for the symmetrical case Eiq = s^q = £15 = £25 = £26 = £35 = £^35 = 
e. The results are plotted in terms of the dimensionless variable (T^ — T)/6T defined as: 



ST T ° \^ 4c, ^ 

Here 6T is the width of the melting transition 

AS + log (^) 

Tm is the melting temperature (neglecting the trimers and higher order structures) for which 
the fraction of rods in the dimer phase F = (2C2) /{Ci + 2C2) = 1/2. 

Tm = ^^7^^-^ (16) 

A5 + log [^^) 



The concentrations for of all the n-mers Cn are 

Ci = 5^Q (17) 

i 

C2 = 5^Q, (18) 

k>j>i 

l>k>j>i 

C5 = 2_^ Cijklm (21) 

m>l>k>j>i 

^6 ^ / ^ Cijklmn \ ) 

n>m>l>k>j>i 

where each index runs over the set {1,2,3,4,5,6}. For the summations it is understood 

that the set of indices must form a connected diagram. For example, the term C1234 does 

not appear in the sum for C4 since the vertex architecture (see Fig. [3]) stipulates that this 

diagram represents two disconnected dimers, C12 and C34. The mass fraction of the n-mers 

M„ is then defined as 

Mn = ^^. (23) 

EkCk 

k=l 

At low temperatures the dominant structure is the 6-mer, which is the fully assembled 
cage except for the binding of rod R4 at vertex D of the cage. With this information at 
hand, we are now in a position to determine the melting profile for the full system, DNA 
linkers together with the DNA grafted nanoparticles. 

IV. DNA-CAGED PARTICLES 

In this section we determine the concentration of nanoparticles decorated with DNA 
structures. We can determine the concentration of particles decorated with DNA structures 
by applying the same rules for the chemical potential as before. For example, consider 
decorating a free particle with concentration Cp with the monomer Ci. We have 

The effective binding energy gi has two contributions. The first comes from the hybridiza- 
tion free energy gi of the DNA arms on the particle hybridizing with the ssDNA ends of rod 

10 
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FIG. 5: (Color online). The mass fraction for the partially assembled cages which can form in 
the absence of the nanoparticle. Mi, M2, M3, M4, M5, and Mq are the mass fractions for the 
monomers, dimers, trimers, tetramers, 5-mers, and 6-mer respectively. In the plot c^°* = InM. 

Ri- As before these hybridization free energies can be decomposed as a sum of contributions 
from joining two arms at a vertex. Analogous to the definition of e, we let g denote the 
average hybridization free energy for joining two rods at the vertex, one of which came from 
the DNA grafted on the nanoparticle surface. 

The second contribution is an entropic contribution associated with localizing the DNA 
structure on the surface of the particle. Since there are Narms DNA strands grafted onto the 
particle surface, there is a combinatorial factor associated with the number of ways to make 
the first connection between the particle and the DNA structure. Let a = N arms / i^'^'^'^) 
be the average areal grafting density of DNA on the nanoparticle surface for a particle of 
radius r, and h ~ Inm a localization length. The entropic contribution can be estimated 
in terms of the concentration if) = a/h which relates g to g in the following manner. 

g = {25r, + 25r, + 5R,)g - T log 



N„. 



^ 



Co 



(25) 



The factor Sr = 1 ii rod Rj is present in the structure, and 6r. =0 otherwise. Here A^ is 
the number of vertices of the cage to which the nanoparticle is bound. In our example case 
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we have N = 2. Putting everything together we have 

Cpi = -^ exp ( ^ j (26) 

gi = 2g-Tlog(Narms^j (27) 

The same general procedure can be apphed to decorating the particles with all of the 
DNA structures considered in the previous section, making sure to take into account the 
vertex architecture. For example, we cannot decorate a particle with the dimer C25 since at 
each of the vertices the DNA arms which come from the particle cannot directly hybridize 
to the arms which come from the dimer. The concentration of particles decorated with 



n-mers is c 



(n) 
■P ■ 



cf = c, ^ (28) 

4^^ = ^.E-^-pf#') (29) 



i 



-?' =c,5^^exp(^) (30) 



cf =c, J^^expMf 
k>j>i ^° ^ ^ 



, J^^expMf (31) 



c^^' =c, Y: ^exp^) (32) 



l>k>j>i 

C 

Co "^ V T 



cf = Cp Yl ^^^^exp(^l^) (33) 



m>l>k>j>i 

4 = Cp 2^ -^^ — exp( ^ 1 (34) 

n>m>l>k>j>i ^ ^ 

If Cp°* is the total initial particle concentration, we can write the equation for particle con- 
servation in the following form: 

cr = E4"^+o(c;) (35) 

n=0 

This equation is then solved to determine the concentration of free particles Cp and hence 
the concentration for particles decorated with DNA structures. The mass fraction m„ for 
particles decorated with n-mers is 



in) 

rrin = -^ • (36) 



Ecf 



-P 
k=0 
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FIG. 6: (Color online). The mass fraction m„ for nanoparticles decorated with n-mers. Note that 
the fully assembled tetrahedral cage surrounding the particle, rriQ, is the dominant equilibrium 
structure for low temperatures. For comparison the mass fraction of the cage in the absence of 
the particles Mq is also plotted. In the plot g = e. 



The results for the mass fraction are plotted (see Fig. [6]) for the case g = e. For low 
temperatures the dominant structure in solution is our desired nanoblock, a DNA-caged 
particle. 

It is of crucial interest for the experimental feasibility of the proposal that the tetrahe- 
drally caged particle is the dominant structure close to room temperature. We can see from 
Fig. [6] that the caged particle is the dominant equilibrium structure for (T^ — T)/6T < —5. 
This in turn determines the standard enthalpy AH and the standard entropy AS for the 
hybridization between two DNA arms at the vertex. We find that AH ~ —100 kcal/mol 
and AS ~ —270 cal K~^ mol~^. For DNA rods with concentration c*°* = InM in a 0.2M 
NaCl solution this gives T^ — 35°C and T ~ 25°C. This information can be used to 
determine the number of DNA bases in each arm of Fig. [2] (i.e. the length of the ssDNA 
ends on the rods). Using the average nearest neighbour parameters of 20|, the DNA arms 
should be 28 base pairs long. Hence each of the dsDNA arms of the 4 arm DNA junction 
is 14 base pairs long. 
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The essential result is that there is an experimentally accessible regime for which the 
dominant equilibrium structure is the DNA-caged particle. In this regime, (T^ — T)/6T < 
—5, we do not expect the assembly of the DNA-caged particle to be kinetically limited. The 
reason is that in this regime, even in the absence of the particle the dominant structure is 
the 6-mer, which is then trivially converted into the DNA-caged particle. 

V. DNA-PARTICLE PARASITES 

In this section we pause to discuss some other DNA-particle structures (parasites) which 
could potentially decrease the overall yield of our nanoblock (see Fig. [7]). One competing 
structure is the particle attached to the outside of a DNA cage. In this case the particle can 
bind to at most three of the tetrahedron vertices. As a result, the equilibrium yield of this 
structure will be proportional to the yield of the DNA-caged particle, but suppressed by a 
factor of the Boltzmann weight for the missing connection exp (^), and thus negligible for 
low temperature. 

We should also consider the stability of our nanoblock with respect to two particle struc- 
tures, i.e. including terms of order c^ in Eq. [35l For example one can consider a particle 
pair attached to the same DNA cage. Building such a structure does not necessarily cost 
binding energy with respect to the DNA-caged particle. However, there is still a cost is 
associated with the loss in translational entropy of the free particle T log ( ^ ) • As a result 
these particle pairs will be relatively rare and should not decrease the overall yield of the 
DNA-caged particle. 

In principle, a particle may have several DNA cages assembled around it. As it turns out, 
these are the most problematic of the competing structures. If we denote the concentration 
of the cage C123456 = Ccage and the free energy for the particle binding to the cage 5^123456 = 
dcage then Concentration of particles with m cages Cp is 

Cf ) = c, (^) ^ exp ^ -^i9cage+^{m-l)a] ^^ ^^^^ 

Here a is an energetic parameter which characterizes the interaction between two cages 
attached to the same particle. Since the dsDNA rods (the cage edges) are charged, this 
interaction is presumably dominated by the electrostatic repulsion o: 



the rods. Within the 
Debye-Huckel approximation this problem has been treated ((2l| , |22|] , 23| ) . The electrostatic 
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FIG. 7: (Color online). Some of the undesired structures which can form. A) One particle attached 
to the outside of a cage. B) Two particles attached to the same cage. C) One particle surrounded 
by two cages. 

energy E{R, 6) of two rods separated by a minimum center to center distance R which make 
an angle 6 when viewed along R is 

E{R,e) f2n\B\exp{-KR) 



T \ kP I sin^ ^^^^ 

2 

Here the rods have the same effective linear charge density z/ = |, A^ = ^ is the Bjerrum 
length, and k,~^ = l/i/lTrXgn is the Debye screening length for monovalent counterions of 
concentration n. Assuming that the electrostatic energy for the cage-cage interactions can 
be expressed in terms of a pairwise sum of contributions from rod-rod interactions we have 



Nr 



a 



{E{R,9)). 



(39) 



Here the angular brackets denote the average, and N^ is the number of rod-rod contacts 
between two cages. An energetically favorable orientation of the cages has N^ = 6. To 
suppress the formation of particles with two cages, we require the following ratio to be small 



Op 



-^cage 



a 



(1) 



exp 



-[gcage + a] , ^^ 



(40) 
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The electrostatic energy can be quite significant. For perpendicular orientations of the 



dsDNA rods reference 2^ reports a contact potential E (rf, ^) ~ SOT (here d ~ 2.4nm is 
the dsDNA diameter) in n = 0.005M NaCl. At fixed temperature and salt concentration, 
Eq. Ho] imposes a condition on the DNA linker concentration which must be met in order 
to suppress the assembly of more than one cage around the particle. If the DNA linker 
concentration is not too high, and the salt concentration fairly low, the assembly of more 
than one cage around the particle can be prevented. 

VI. HIERARCHICAL SELF-ASSEMBLY 

The DNA-caged nanoparticles in this paper are interesting nanoblocks in a hierarchical 
self-assembly scheme. Part of their usefulness stems from the fact that interactions between 
nanoblocks are highly anisotropic. Recall that at each vertex of the DNA cage there is a 
unique ssDNA sequence available for hybridization. As a result two nanoblocks can be made 
to interact in a very specific manner by introducing another set of vertex- vertex DNA linkers. 
Moreover, the number of these vertices is explicitly determined by the cage architecture, 
which translates into a well defined "valence" for the interactions between nanoblocks. 

Here we discuss a particular hierarchical self-assembly proposal where DNA-caged par- 
ticles are the natural building blocks. The basic proposal of the Voronoi scheme is the 
following. Any target structure not necessarily crystalline can be represented as a discrete 
set of points, i.e. the location of particles in the structure. With this set, one performs the 
Voronoi decomposition |24]. The Voronoi cells for the structure can be used to design the 
cages surrounding the particles in the following manner. For any given particle, place a 
vertex at the midpoint between that particle and each of its Voronoi neighbors. In this way 
we map the target structure onto the set of caged particles, with certain pairs of vertices 
to be connected. These vertices can be connected by introducing a set of DNA linkers 
which perform the vertex- vertex binding. By construction the target structure must be the 
ground state of the system. What remains to be seen is whether or not the target structure 
is favored kinetically. A task for future research is a detailed consideration of this question 
using Monte Carlo simulations. 
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FIG. 8: (Color online). A simple example of the Voronoi scheme in two dimensions. 



VII. CONCLUSIONS 



In this paper we discussed a proposal to self- assemble DNA-caged particles. The basic 
components are several types of rod-like dsDNA linkers with ssDNA ends, and nanoparticles 
grafted with ssDNA. By designing the ssDNA sequences appropriately, the dsDNA rods 
self-assemble into a cage surrounding the particle. The edges of the cage are dsDNA, 
and the vertices are multi arm DNA junctions. A particular implementation of this idea 
was discussed for the self-assembly of tetrahedrally caged particles. We calculated the 
equilibrium yield of the DNA-caged particles and discussed their stability with respect to 
alternative structures. At low temperature, the nanoparticle surrounded by one cage is the 
dominant equilibrium structure. Although the calculations were performed for a tetrahedral 
cage geometry, the ideas are generally applicable to many types of polyhedral cages. Each 
vertex of the cage with degree V can be constructed from a l^ + 1 arm DNA branched 
junction. Such junctions have been constructed with up to 12 arms, which leads open the 
possibility of much more complicated cages ((25|,|26|). A natural next step would be to 
consider a particle inside a DNA cube, since the same vertex architecture proposed in this 
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paper would apply. 

We concluded by discussing the usefulness of DNA-caged particles in a hierarchical self- 
assembly proposal. The Voronoi scheme maps the problem of self-assembling a particular 
target structure onto a set of caged particles for which the target structure is the ground 
state. The experimental realization of self-assembled DNA-caged particles would represent 
an important step towards realizing the technological potential of DNA. 

Thanks to A. Klopper for help with the figures. This work was supported by the ACS 
Petroleum Research Fund (Grant PRF No. 44181-AClO). 
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VIII. APPENDIX: OVERLAP DENSITY Ce// 

In this appendix we calculate the overlap density Ce// used to determine the effective 
hybridization free energy for the DNA structures in Fig. HI It is helpful to consider the 
associated problem of determining the probability distribution for the end vector of a freely- 
jointed chain made up of A^ linkers, each of length L. The probability distribution for the 
chain composed of one linker is simply. 



6{R - L) 



P.(R) = '^^^^^ (41) 



Therefore for a chain composed of A^ such linkers we have: 



Pn 



;R) = n / d%p,{r,)6' (y^ r, - R j (42) 
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The inverse Fourier transform of the probabihty distribution has a particularly simple form 

sm{kL)^ 



PN{k) = / d^Rexpitk ■ R)pjv(R.) 



kL 



PnCR) = {2n)~^ / d3kexp(-ik ■ R)p7v(k) 
Working in spherical coordinates, performing the angular integration yields 

_3 r,,fsm{kL)Y 



pjv(R) = (27r) 



kL J 



dk / dilexp{—ikR cost 



(43) 
(44) 

(45) 



2 /•°°. . ,. _ /sin(A;L)\ 
kL ) 



N 



(27r)~^- / ksm{kR) 

R Jo 



dk 



Here Jq{z) = ^^^^ is the spherical Bessel function 281] of order 0. 

We are now in a position to determine the overlap density Ce// for the triangle structure. 
Of particular interest for this calculation is ( 29|, 301]): 

e{2L - R) 



P2(R) 



StiL'^R 



(46) 



Here 6{x) is the Heaviside step function. The overlap density for the triangle is then 
calculated as 



-e// 



d^K / dhp2iK)pi{r)6%R-r) 



(47) 



1 



87rL2 47rL2 



R 



327r2L4 
1 



dh 



e{2L - r) 



5(r - L) 



SttLS 



There are two remaining overlap densities which need to be calculated. One of them is 
for forming an equilateral parallelogram, i.e. a diamond. In Fig. H] this structure is the 
third structure in the tetramer row. Forming this structure can be viewed as joining the 
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ends of two chains, each of which is composed of two hnks. 



Ceff = fd^Rf dhp2{R)p2ir)5%R - r) (48) 



1 V /".So f,.0i^L-R)e{2L-r) 



— — I / rf'R / d'^r"- — ^-^^- -6%R 



.e{2L-r) 



d-'r 



1 



87rL3 
Alternatively, one can calculate the overlap density for the diamond structure as 



"e// 



ld'Rldhp,{K)p,{r)6'{R-r) = ^ (49) 



This is simply viewing the diamond as joining a chain of three links with a chain of one link. 
The results are the same, as they must be. 

The remaining overlap density to be calculated is associated with making the last con- 
nection in the tetrahedron. Assume that all the connections have been made, except for 
the connection between rod R4 and vertex A. Taking vertex D as the origin, the position of 
vertex A is the end vector of a two segment chain, with each segment having length / = ^^y^. 
With one end of rod -R4 fixed at vertex D, the other end must connect to vertex A. 

Ceff = f d^R f dhp2(R)piir)6%R-r) (50) 



d'R / d-'r-—- — '-5{r -L)5' (R-r) 



SttP 47rL2 J J R 

' d'^r— -6{r - L) 



327r2/2L2 
1 



SnPL 
1 

67rL3 
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